Adaptive Zeno measurement for qubit state preparation 
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This paper reconsiders the method of adaptive measurement for qubit state preparation developed 
by Jacobs and shows an alternative scheme that works even under unknown unitary evolution of the 
state. The key idea is the use of quantum Zeno effect; i.e., the measurement is adaptively changed so 
that the state is dynamically frozen in one of the eigenstates of the measured observable while that 
eigenstate approaches to the target state. Numerical simulations demonstrate that this adaptive 
Zeno measurement scheme offers faster and more stable convergence of the state to the target as 
well as higher robustness against the uncertainty of the unitary evolution, than the case of Jacobs' 
scheme. 
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I. INTRODUCTION 

Repeated measurement of an observable that is appro- 
priately changed according to the pre-measurement out- 
comes, i.e., the adaptive measurement, has great poten- 
tial for various purposes in quantum information sciences. 
The first demonstration has come out in the application 
to quantum phase estimation [3-[3]. Another application 
of adaptive measurement is deterministic state prepara- 
tion [5|-|7|, which all consider the case of a single qubit. 
A striking feature of this scheme is that a desired time 
evolution of the state is brought only by measurement 
back-action, and there is no need to introduce any exter- 
nal force for controlling the state. 

Let us especially focus on the method developed by 
Jacobs 0. This employs the schematic of a continuous- 
time measurement; in this case the probabilistic change 
of a qubit state p is described by the following stochastic 
master equation (SME) 0,0: 

dp t = - k t [a t , [a t ,pt]]dt 

+ y/Zkt{& t pt + Pt&t ~ 2Tr {a t pt)pt)dW t , (1) 



where dW t is the standard Wiener process satisfying the 
Ito rule dWt — dt. In the above equation, b t and k t 
represent the measured observable and the measurement 
strength, respectively. Adaptive measurement means 
that we can change at and kt continuously in time, as 
functions of the state pt, so that pt will converge to a 
target state. In Jacobs' scheme, the state is assumed to 
be pure (this actually holds if the initial state is pure) 
and is thus of the form 

p t = \^ t )(vjtl |^>=cos(<5 t /2)|0)+sin(5 t /2)|l}, (2) 

where |0) = (1, 0) T and |1) = (0, 1) T are the target state 
and the initial state, respectively. Then the observable 



and the measurement strength are updated according to 
the following laws: 



a x cos(St) - a z sin((5 t ), fc t J 



kS?, 



(3) 



where k is a positive constant. This means that the mea- 
sured observable a\ is changed so that its eigenstates 
are always perpendicular to the current state \ipt) in the 
Bloch sphere representation, as shown in Fig. 1. The 
measurement strength kj is also adaptively changed, and 
it decreases proportionally to the distance between the 
current and target states; hence the measurement be- 
comes very weak when the state approaches to the tar- 
get. In fact, with this adaptive measurement law <j3j) , the 
time evolution of St £ [0, tt] is given by 



dS t = 



^KStdWt, (4) 
that S t converges to 



and it was numerically shown in 
zero as t — > oo almost surely. 

Here we come up with the question about how much 
the above adaptive measurement scheme is robust against 
certain disturbance acting on the system. Actually such 
robustness against a specific decohering effect was eval- 
uated in Q. Another specific but important disturbance 
is an unknown unitary evolution of the qubit state, in 
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FIG. 1: (Color online) The measurement axes of Jacobs' 
adaptive measurement scheme (left) and the adaptive Zeno 
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which case the driving term — i[H, pt]dt is added to the 
right hand side of Eq. ([!}. For example, if we take a 
two-level atom continuously observed using the Faraday 
rotation technique (To| to realize the qubit system sub- 
jected to Eq. ([T]), such a disturbing Hamiltonian may 
appear and take the form H — Aa y with A unknown de- 
tuning between the atomic transition frequency and the 
laser frequency. Note that in [5( this kind of unknown 
disturbance was not discussed. 

In this paper, we reconsider the same control problem 
discussed above, yet with additional care about the influ- 
ence of an unknown unitary evolution, and then propose 
a new adaptive measurement scheme that has clear ro- 
bustness property against that disturbance. A novel fea- 
ture of this scheme is in the use of quantum Zeno effect 
[HHIl], which states that a repeated (or continuous in 
the present case) measurement freezes a quantum system 
in a particular state, under a unitary evolution. There- 
fore, it is fully expected that our adaptive measurement 
scheme having the Zeno effect can keep the system stable 
around an arbitrary target state even under influence of 
an unknown Hamiltonian added to the system. It will be 
actually shown by numerical simulations that, when the 
state approaches to the target, the quantum Zeno effect 
certainly emerges, suppresses the disturbance, and sta- 
bilizes the state; as a result, our adaptive measurement 
scheme offers faster and more stable convergence of the 
state to the target as well as higher robustness against 
the uncertainty of the unitary evolution, than the case of 
Jacobs' scheme. 



II. ADAPTIVE ZENO MEASUREMENT 

When we measure a system observable enough fre- 
quently, then the state is frozen in one of the eigenstates 
of that observable even under a certain unitary evolution. 
This is called the Zeno effect. Hence, if this measured 
observable is changed adaptively so that the eigenstate 
dragging the current state approaches to a desired target 
state, it is expected that the state will finally be stabi- 
lized at that target state. Based on this idea, the fol- 
lowing observable is taken as a measured observable in 



Eq. Q: 



sin(aiJt) + a z cos(ai5t), 



(5) 



where a € [0, 1] is the parameter whose meaning is ex- 
plained below. The eigenstates of the observable af are 



, _ / cos(a£ t /2) 
' z) - I sin(aJ t /2) 



-z = 



sin(a£t/2) 
- cos(a(5 4 /2) 



which satisfy &f\+ z ) = |+z), ^ t z |- z ) = -|- z >. In 
the Bloch sphere representation, |+z) divides the an- 
gle between the target state |0) and the current state 
\ij} t ) = (cos((S t /2), sin((5 t /2)) T into a : (1-a); see Fig. [Q 
The transition probabilities of the state jumping to these 
eigenstates are p + = |(+z|V>t)| 2 = (1 + cos(/36 t /2))/2 and 



P- = l(-z|^ t )| 2 = (1 - cos(/35 t /2))/2, where /? := 1 - a, 
and they satisfy p+ > p_ . This means that the state 
tends to move towards the direction of |+z) with prob- 
ability p + . At the same time, since in this case St de- 
creases, the eigenstate \+z) moves towards the target |0). 
In particular, when the state reaches the target |0), or 
equivalently St — > +0, the eigenstate \+z) also converges 
to the target and the transition probability p+ takes the 
value 1; that is, the Zeno effect emerges, and the state is 
frozen in the target as expected. For this reason, we call 
this method the adaptive Zeno measurement scheme or 
simply the Zeno measurement scheme in what follows. 

Next, the measurement strength is simply set to a con- 
stant value kf = fc z , unlike the case of Jacobs' scheme. 
The reason will be explained in the remark (i) given in 
the end of this section, but we here point out that this has 
a clear merit from a practical viewpoint. In fact, chang- 
ing the measurement strength in addition to changing 
the measured observable definitely costs more expensive 
compared to the case where only the latter is required. 
In this sense, our scheme is suited to experimental im- 
plementation. 

Regarding the disturbing Hamiltonian, as discussed in 
Introduction, we take H — Aa y , where A is an unknown 
constant. Note that with this Hamiltonian the state ro- 
tates around the y axis in the Bloch sphere. 

Consequently, the dynamical evolution of the state 
under the Zeno measurement setup and the disturbing 
Hamiltonian introduced above is given by 

dp t = - iA[& y ,pt]dt- k z [af, [cr?,p t ]}dt 

+ V2kZ(&fpt + ptSrf - 2Tr (afp t ))dW t . (6) 

The initial state is now on the x-z plane, hence we have 
the dynamics of St as follows: 

dS t = 2Adt - 2fc z sin(2/3(5 t )di + V8k% sm(P5 t )dW t . (7) 

The linear approximated equation, which is valid when 
S t ~ 0, is 



dS t = 2Adt - 4:k z f3S t dt + V8k z /35 t dWt. 



(8) 



On the other hand, the dynamics of the same St but with 
Jacobs' scheme is given by 



dS t = 2Adt 



'SKStdWt. 



(9) 



We now compare the dynamical stabilities of the above 
two dynamics ([SJ and (JSJ , in an intuitive way. To do this 
let us ignore the boundary condition St € [0, tt] and also 
set A — 0; then Eqs. and © can be both analytically 
solved as follows 1141: 



(Z) 5 t =exp (~4k z 'l3-Ak z l3 2 )t + V8k z (3Wt S a , (10) 



(J) St 



exp 



Ant 



<&KWt 



So, 



(11) 



where Wt = J dW s . A fair comparison of the dynamical 
stabilities of these equations can be performed by taking 
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the same values of diffusion coefficients; i.e. 



V8k, which leads to k /3 2 = K. In this case, Eq. (JTUJ) 
becomes 



(Z) 5 t = exp 



4 1 



Kt 



inWt 



S Q . (12) 



Equations (jlip and (fT2|) imply that, when the two adap- 
tive measurement schemes introduce the back-action 
noise of the same magnitude, the Zeno measurement 
scheme offers the faster convergence of the state to the 
target of order e -4 "*/^; this term indicates no more than 
the emergence of the quantum Zeno effect. Of course, 
the above discussion makes sense only when the state is 
around the target. Thus we need numerical simulations 
to actually compare the stabilities of the two schemes, 
which will be discussed in the next section. 

Before closing this section, we provide two remarks. 

(i) Jacobs' scheme requires the adaptive tuning of the 



measurement strength (i.e., k 



J _ 



c<5 t 2 in Eq. ©) for 



the dynamics of St to have the state-dependent diffu- 
sion term. As seen from the explicit solution (fTTj) . such 
state-dependence is indeed necessary to generate dynam- 
ical stability of S t . Hence, it should be maintained that, 
with the Zeno measurement scheme, the diffusion term 
of Eq. ([7]) depends on the state even with the fixed mea- 
surement strength (i.e., kf = k z ). 

(ii) The time evolution of the fidelity between the cur- 
rent state pt and the target state |0) is given by 



d(Q\p t \Q) 



A sin St + k z ^ cos(2a<5t — St) — cos St 

3 + cos 5 t , . . , 
cos(ad t — S t ) + cos(ad t/ 



dt 

dW t . 

(13) 



Hence, the optimum value of a that maximizes the deter- 
ministic change per unit time of the fidelity is given by 
a = 1/2. This value is actually taken in the simulations 
shown later. 



III. NUMERICAL SIMULATION 
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FIG. 2: (Color online) Time evolutions of (a) the mean and 
(b) the standard deviation of St when A = 0. The solid blue 
and the dashed red lines correspond to the Zeno measurement 
and Jacobs' scheme, respectively. 
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FIG. 3: (Color online) Time evolutions of (a, c) the mean 
and (b, d) the standard deviation of St for several values of A. 
Figures (a) and (b) show the cases of the Zeno measurement, 
while Figs, (c) and (d) for the Jacobs' case. The green, 
red, water blue, purple, and blue lines correspond to A = 
1,1/10, 1/10 2 ,1/10 3 , and A = 0, respectively. 



In this section, we numerically compare the perfor- 
mance of the Zeno measurement scheme to Jacobs' one, 
under the assumption that A is known; this is a crucial 
assumption because we are then able to update St and 
as a result the adaptive measurement law (J3J) (or Eq. 
([3])) exactly by recursively solving Eq. ([7]) (or Eq. © 
for Jacobs' case). The parameters are taken as n = 1 
and k z = 4. Also we set a = 1 — j3 = 1/2, based on 
the fact mentioned in the remark (ii) in Sec. II. These 
parameters satisfy the condition for the fair comparison, 
k z (3 2 = k, which was given above Eq. (|12l) . The dis- 
turbance strength A takes several values. The initial 
condition is S — n, as defined below Eq. ([2]). 

Fig. [2] shows the time evolutions of the mean and the 
standard deviation of St, when A = 0. The plots are 



obtained by averaging 10 6 sample paths. In Fig. [5] (a) 
the solid blue and the dashed red lines correspond to 
the cases of the Zeno measurement and Jacobs' schemes, 
respectively. Both lines can be fitted by a simple func- 
tion F(t) — Cie~ C2t ; the decay rates are obtained as 
C-2 = 10.03k for the former case and C2 = 1.357k for the 
latter case (note k = 1 now). That is, the convergence 
speed in the Zeno measurement case is about 7.4 times 
faster in average than that in Jacobs' case. Moreover, 
it can be seen in Fig. [5] (b) that the state evolves with 
smaller fluctuation, or in other words the time-evolution 
of the state is more stable, when using the Zeno measure- 
ment scheme. Note that the Zeno measurement scheme 
was originally expected to offer faster and more stable 
convergence of the state to the target when A ^ 0, but 
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it is indeed the case in the ideal case A = as well. 

Next we examine several values of A and show in Fig. [3] 
the corresponding means and standard deviations of St. 
Figure[3](a, c) shows that the mean (St) does not converge 
to zero when A ^ but rather takes a constant error in 
the long time limit; in particular, we numerically find 
that this constant error is given by e = 0.2544Ak for the 
Zeno measurement case (a) while in Jacobs' case (c) there 
is no such a simple rule. Clearly, the Zeno measurement 
scheme moves St towards the target 5 = faster with 
smaller constant error in average than Jacobs' case, for 
every values of A. Similar plots are observed in Fig.[3](b, 
d); that is, the state approaches to the target with smaller 
fluctuation, when the system is governed by the Zeno 
measurement scheme. 



IV. ROBUSTNESS OF THE ADAPTIVE ZENO 
MEASUREMENT 

Here we study the case where the external force A 
is unknown. To make the situation clear, let us again 
consider the general SME ([l} that is additionally driven 
by an unknown Hamiltonian H: 

dp t = - i[H, p t ]dt - k t [a t , [a t , Pt]] dt 

+ yfZk t [atpt + Pta t ~ 2Tr (a t p t )p t ]dW t . (14) 

This true state pt cannot be precisely updated, due to 
the uncertainty of H. Therefore, we need to devise an 
updating law of a nominal state, say p' t , only using the 
measurement result yt that is subjected to the output 
equation 

dy t = Tr (a t Pt )dt + dW t . (15) 

Note this is driven by the same dWt as that in Eq. (fl"4]) ; 
dWt is called the innovation in the framework of quan- 
tum filtering theory [15l 1 1 61 ] . As an updating law of p' t , 
we particularly use Eqs. (Ti"4|) and (fl~5|) with p t and H 
replaced by p' t and H', respectively: 

dp' t = - i[H', p' t ]dt - k t [at, [a t ,Pt]\dt 

+ yfZk t [a t p't + P'A - 2Tr (a t p' t )p' t ] 

x[dyt-Tr(a t p' t )dt\. (16) 

Note again that y t is the measurement result and is thus 
known. Hence we can recursively calculate the nominal 
state p' t , although it should differ from the true state p t . 
In particular, in the adaptive measurement setup, the 
observable a t and the strength k t are changed in time as 
functions of p' t . 

For the specific problem under consideration, let H = 
Aaz and H' = A'az be the true and nominal Hamilto- 
nians, respectively; A is an unknown constant while A' 
is a known nominal constant. Also we define St and S' t , 
corresponding to the true and the nominal states. In the 
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FIG. 4: Sample paths of the true variable St (solid blue lines) 
and the nominal one S' t (dashed red lines) for (a) the Zeno 
measurement scheme and (b) Jacobs' scheme, where A = 
1/10 2 , A' = 0. 

Zeno measurement scheme, these variables are driven by 
the following equations: 

dS t = [2A - 2k z sm(28 t - 2aS' t )]dt 

+ V8k^sia(St-aS' t )dWt, (17) 
dS' t = [2A' - 2k z sin(2j8$)] dt 

+ V8k? sm(j38't)dWf , (18) 

where 

dWf = [cos(/3S' t ) ~ cos{S' t - aS t )]dt + dW t . 

Note that j3 = 1 — a. The above two equations (fTTj) 
and (I18p take the same form when A = A' and S t = S' t . 
Jacobs' scheme, on the other hand, leads to the following 
equations to update St and 5[: 

dS t = [2A + 2kS[ 2 sin(2<5 t - 2S' t )]dt 

+ V8iz5' t cos(S t - S' t )dW t , (19) 
dS' t = 2A'dt + V$^S' t dW t 3 , (20) 

where 

dWl = sm(S' t -S t )dt + dW t . 

Again, the above two equations (fT9"| and (|2U|) take the 
same form when A = A' and St = S' t . 

In Fig. [5]we show sample paths of St (solid blue line) 
and S't (dashed red line) for both of the adaptive mea- 
surement schemes, with the same parameters taken in 
Sec. Ill, fixed uncertainty A = 1/10 2 , and a typical nom- 
inal value A' = 0. It is apparent from the figures that, 
with the Zeno measurement scheme, the true variable S t 
is stabilized around the target S = though it still has 
a certain error corresponding to the constant external 
force, while in Jacobs' scheme St finally follows a deter- 
ministic time-evolution and goes away from zero; in every 
numerical simulation a similar behavior is observed. 

We here address the mechanism that brings about the 
above drastic difference. First, since the nominal state 
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FIG. 5: Sample paths of the true variable St (solid blue 
lines) and the nominal one 8[ (dashed red lines) for (a, c, e) 
the Zeno measurement scheme and (b, d, f) Jacobs' scheme. 
For all cases, the true external force is set to A = 1/10 2 , while 
several nominal values of A' are examined: (a, b) represent 
the case when A' = 1/10 3 , (c, d) for the case A' = 1/10 2 , 
and (e, f) for the case A' = 1/10. 



is now governed by the dynamics without uncertainty, 
S' t converges to zero, as seen in Sec. III. When S' t — > 0, 
the measured observable becomes at — o z for both the 
adaptive measurement schemes, but the measurement 
strength becomes k( — for Jacobs' case while kf = k z 
for the Zeno measurement case; that is, Jacobs' scheme 
stops to measure the system once S' t = is reached, and 
consequently the true variable St turns out to follow the 
deterministic dynamics dSt/dt — 2A. In contrast, in the 
Zeno measurement scheme, even after S' t takes zero, a z 
is continuously measured, implying that the true state is 
stabilized around the target. 

The above statement holds only for the case of A' = 0; 
hence we take some non-zero values of A' and show in 



Fig. [5] the trajectories of St and S' t . If A' is not zero but 
smaller than the true value A (Figs, (a, b)), we observe 
similar trajectories to the case of A' = 0; that is, in Ja- 
cobs' case the true variable again goes away from zero 
deterministically, due to the same reason stated above. 
In Figs, (c, d) where A = A', the nominal variable S' t 
exactly tracks St, which is the property the nominal up- 
dating laws (US]) and (P2U1) should have. Lastly, if A' is 
bigger than A (Figs, (e, f)), or in other words if we take a 
conservative approach for estimating St, Jacobs' scheme 
can stabilize the state around the target although it has 
a big fluctuation, while the Zeno measurement moves the 
state more close towards zero with much smaller estima- 
tion error. 

Summarizing, in the numerical simulations shown here, 
the Zeno measurement scheme always offers smaller es- 
timation error between St and 5[. That is, it is fairly 
robust against the uncertainty of the Hamiltonian. This 
nice property should be particularly maintained, because 
without robustness the stabilization can fail as seen in 
Figs.H(b) and [5](b). 



V. CONCLUDING REMARKS 

The main features of the Zeno measurement scheme 
are twofold; the first is that the measured observable con- 
tains an eigenstate that eventually converges to the tar- 
get state, and the other is that the measurement strength 
does not become weakened even when the state is close 
to the target, which is the crucial property bringing the 
robustness property. As long as the above two conditions 
are satisfied, it is expected that the Zeno measurement 
scheme can be generalized to the multi-dimensional case. 
However, in many situations any physically available ob- 
servable does not contain an eigenvector that is identical 
to a desired target state, and then the Zeno measure- 
ment scheme cannot be used in order to prepare that 
target state. Rather, as discussed in in such a case, 
a certain weak measurement may work for generating a 
target state. But a critical requirement for this measure- 
ment is that the measurement strength should be weak- 
ened when the state approaches to the target; then, the 
system can become fragile against unknown disturbance, 
as seen in Sec. IV. Exploring an adaptive measurement 
method that overcomes this issue is an interesting future 
work. 

We thank Yuzuru Kato for his contribution in the re- 
mark (ii) in Sec. II. 
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